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Abstract. We present imaging and spectroscopic observations with HST (WFPC2, ACS/HRC and STIS), VLT (FORS2) 
and Keck (HIRES) of the dM8 ultra-cool dwarf DENIS-P J 185950.9-370632, located in the R-CrA region. The presence of 
lithium absorption at 670.8 nm and the strong Ha emission indicate a young age and a sub-stellar mass. Our diffraction-limited 
images resolve a companion at the separation limit of HST/ ACS (~ 0'.'06). The 2.1 mJy flux in the LW2 filter (5.0-8.5yum) 
of the Infrared Space Observatory jOlofsson et'an[l999 ) likely corresponds to an infrared excess, suggesting the presence of 
circumstellar material. Proper motion and photometric measurements, as well as the Ha activity, confirm membership in the 
R-CrA star forming region. If confirmed by further observations, DENIS-P J 185950.9-370632 would be the first accreting 
sub-stellar multiple system observed to date. 

Key words. Stars: individual (DENIS-P J185950.9-370632) - Stars: low mass, brown dwarfs, disk - Binaries: visual - 
Techniques: high angular resolution 

Note: Version with higher resolution figures available on molecular cloud clumps ('Bodenheimei'igggl'Shu et al.'l987l), 
\ftp://ftp. mpe.mpg.de/people/hbouy/publications or do they form more like "planets", in gravitationally unstable 

regions of circumstellar disks jPapaloizou & Terauemll200i ? 
The existence of brown dwarf binaries or the presence of disks 
are essential diagnostics of these formation scenarii. 



1. Introduction: Disks and Binaries at the 
sub-stellar limit 



Over the last few years large numbers of sub-stellar objects Several dozens very low mass multiple systems have been 

have been discovered as a result of the new sky surveys ^^nd both in the field and in open clusters. In the field the bi- 

(2MASS, DENIS, SDSS, see respectively Skrutskie 1997^ nary fraction^Jor^separations betweenSand 20 A.U is -15% 

rEDchteinetalJll99l: making detailed studies of '^ouv et al. 2003; Burgasser et alj |2QQ3 I Closp gt^lj |2QQ3 

their physical properties possible. One of the most debated »Gizis et al., 2003J. A recent study of colour-magnitude dia- 

question is how do brown dwarfs ori ginate? Are they ejected S^ams in the Pleiades suggests that the overall binary fraction in 

stellar embryos t Reipurth & Clarkel |ooi iBateetalJ l20o3 oP^n cluster might be as high as ~ 50+!/,% ( Pinfield et al.l 

iDelgado-Donateet al.. .20031) ? Do they instead form isolated, 12003), although recem imaging surveys (.MjaULeUd .. 2()03; 

like ordinary stars, due to the fragmentation of collapsing ^^id & Mahond|2Q00t | Basri & Martfn| | l999|) found at most 

~12% visual binaries (separations greater than ~7.5 A.U). 

Send ojfprint requests to: H. Bouy Surveys looking for disks have also found large numbers of 
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disks around very lo w mass objects in star forming regions 
and clusters (see e.g IPascucci et al.ll2003t llfesti et al.ll2002t 
Javawardhana et al .'2002': Natta & Te sti"2001: Muz erolle et alJ 
200(1 , and references therein). From near-infrared observa- 
tions, |j^awardhana et al. (2003) estimated initial frequencies 
of disk around brown dwarfs ranging from 33% to 70% for 
the p-Oph, IC 348, Taurus, C ha I, Upper Sco cr-O ri and TW 
Hya star forming regions, and lMuench et alJ ( l200lh measured 
over 50% in the Trapezium Cluster. Although this suggests that 
brown dwarfs originate like stars more than like ejected em- 
bryos (which should have no disk), our understanding of the 
formation processes is far from complete. The truth might well 
be a combination of the different processes that have been pro- 
posed. 

In the course of a large search for field late-M and L dwarfs, 
(Delfosse et al. in prep.) we have selected for follow-up spec- 
trosc opv all objects in 570 squares degrees of the DENIS sur- 
vey dEpchtein et al.lfl997l) at high galactic latitude that have 
I - J > 3.0 mag. A few fields at intermediate galactic lati- 
tudes (15°to 20°) overlap some nearby star forming regions. 
In the Corona Australis (CrA) molecular cloud complex, we 
selected one red object, with colours typical of very late M- 
dwarfs, DENIS-P J 185 950.9-370632. The CrA region (see e.g 
Neuhiiuser et al. 2000) is o ne of the nearest re gions (~130 pc; 
Marraco & Rvd^rea.1 981.: .Casev et al.lll 998l) where interme- 
diate and low-mas s stars form. Its estimated age ranges be- 
tween 0.5-10 Myr (Knacke et al.'l973; Preibis ch & ZinneckeJ 
1 19991 iNeuhauseret al. .2000) . The complex contains a dark 
cloud close to the R-CrA star, reflection nebulae, and sev- 
eral embedded infrared (IR ) sources , some of wh ich are IR 
class I objects jAdams et alJl 1987; Andre & Montmerl&1994 ). 
A numbe r of T Tauri stars ar e associated with the CrA 
cloud (see N euhauser et al 12000') and the faintest of those are 
brown dwarfs candidates or objects in the transition region be- 
tween low mass stars and br own dwarfs IWilking et al. 1997; 
Ipernandez & Comer6rl 20oT) with spectral types around M6. 

In this paper, we present the results of optical and near- 
infrared imaging and spectroscopy of DENIS-P J185950.9- 
370632. In section |2l we will describe the observations, then 
in section|3]we will give the results of these observations, and 
in section|4]we will analyse these results in the context of re- 
cent models of formation and evolution. 

2. Observations 

2.1. DENIS observations 

The DENIS observations are carried out on the ESO 1 m tele- 
scope at La Silla. Dichroic beam splitters separate three chan- 
nels, and a focal reducing optics provide scales of 3'.'0 per pixel 
on the 256x256 NICMOS3 arrays used for the two infrai-ed 
channels, and I'.'O per pixel on the 1024x1024 Tektronix CCD 
detector of the I channel. The image da ta were processed w ith 
the standard DENIS software pipeline ("Borsenberger" 1997") at 
the Paris Data Analysis Center (PDAC). Source extraction and 
photometry ar e performed a t PDAC, using a space-varying ker- 
nel algorithm ('Alard'"2000'). When searching very rare objects 
in a large data base, the first challenge is artefact rejection. We 



use a set of morphological parameters, based on the correlation 
between PSF model and the object profile and on the consis- 
tency between several aperture magnitudes (Delfosse et al. in 
prep.). 

DENIS-P J185950.9-370632, was selected as a candidate 
very low mass star or and brown dwarf due to its red colors 
(I - J - 3.0 + 0.1 mag). It is located in the R-CrA com- 
plex between the TY-CrA and e-CrA stars (see Figure 
and very close to the core of the ve rv dense molecular cloud 
(Ay ~ 45 mag; WiUcing et al. '1992'). Locally the absorption is 
low, however (see Section ITU . Figure |5] shows finding charts 
in I and J. 

The apparent magnitude of DENIS-P J 185950.9-370632 is 
consistent with a ~5 Myr late-M dwarf at the distance of the 
CrA complex (see section which makes it an excellent 
brown dwarf candidate. 

In addition to the DENIS survey, DENIS-P J 185950.9- 
370632 has been detected in several sky surveys, such as the 
USNO-B 1 .0 catalogue (iMonet et alT2003l , where it is reported 
as USNO-Bl.O 0528-0926219), the GSC2.21 catalog (reported 
as GSC2 S33202002902), and the 2MASS survey, where it is 
reported as 2MASSW18595094-3706313. Tableffland Figure 
|3]summarize these astrometric and photometric measurements. 

2.2. Observation Summary 

Table|2]gives an overview of all the observations we have been 
conducting or retrieved from the archives. We obtained high 
angular resolution imaging using the Hubble Space Telescope 
(HST). Spectroscopy was obtained at both high (Keck-HIRES) 
and low (VLT-F0RS2 and HST-STIS) spectral resolution. 
Finally, we retrieved ISO archival data in which DENIS- 
P J 185950.9-370632 was detected. 

2.2.1. Imaging 

We observed DENIS-P J 185950.9-370632 twice at high angu- 
lar resolution usin g the HST, fi rst with the Planetary Camera 
(PC) of WFPC2 jBaggett et alJ I2OO2I) an d second with th e 
High Resolution Channel (HRC) of ACS d Mack et aljEool . 
The WFPC2 data were taken in SNAPSHOT mode (program 
GO8720, P.I Brandner). We took single exposure for each fil- 
ter (F675W and F814W), to minimise overheads and increase 
exposure times. This prevent automatic removal of cosmic ray 
hits, but we were fortunate enough that no cosmic ray fell close 
enough to the target that it would affect the analysis. The ACS 
data (program G09451, P.I Brandner) were obtained on 2002 
September 24th in CR-SPLIT mode and with a four points 
dithering pattern in each of the F625W, F775W, and F850LP 
filters. This allows correction for cosmic ray events and bad 
pixels. The WFPC2 and ACS data together provide accurate 
photometry in five different optical filters and positions at two 
epochs. The target-acquisition frame for the STIS spectroscopy 
gives a third-epoch position. 

DENIS-P J 185950.9-370632 is obviously elongated on the 
WFPC2, ACS and STIS images (Figure gj. PSF subtraction 
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brings up the companion more clearly on the WFPC2 and ACS 
images (Figure |5}. 

We retrieved ISO archival data of the R-CrA region 
obtained on 1996 April 20th (TDT 15500328, CAMOl, 
[oiofsson 'etaDll999l) with ISOCAMl whic h include DENIS- 
P J185950.9-370632. IOlofsson et alJ ( fl999l) report a 2.1 mJy + 
0.4 mJy detection in the LW2 filter (5.0-8.5yum), at a position 
that is 4'.'0 away from our coordinates for DENIS-P J 185950.9- 
37063 2. This is well within th e pointing uncertainties of ISO 
(~6'.'0, iBlommaertetaP bOOll). and we identify the infrared 
source, ISO-CrA 63, with DENIS-P J 185950.9-370632. The 
source is not detected in the LW3 filter (12-18/vm), and we 
measured a 3 cr upper limit of 1 . 1 mJy. No other optical/near-IR 
source is present in the ~6" ISO error box. 



slit oriented along the axis of the binary. Unfortunately, the 
small separation (~0'.'060), the relatively small flux ratio (~0.3 
in R and I, Section and the relatively low S/N prevent us 
from properly resolving the two components (the pixel scale 
of STIS spectra is ~0'.'050). We processed the integrated spec- 
trum using the recommended STSDAS tools in IRAF and the 
best calibration files provided by the STScI archive. 



3. Analysis 

3.1. Spectral Type and extinction 



2.2.2. High Resolution Spectroscopy 

DENIS-P J185950.9-370632 was observed on 30 May 2000 
with the High Resolution Echelle Spectrometer (HIRES; 
IVogt et al . 1994) on the Keck I telescope. A slit width of r.'15 
and two-pixel binning in the spectral direction gave a resolving 
power of R=33 000. The exposure time was 1 200 s and the air- 
mass 1.83. Fifteen echelle orders were recorded on the detec- 
tor, covering the wavelength range from 667. 1 nm to 895.0 nm, 
with gaps between the orders. The data were reduced following 
a standard procedure in the IRAF' environment, including bias 
subtraction, flat-field correction, aperture extraction and wave- 
length calibration using a Th-Ar lamp spectrum. 

2.2.3. Low Resolution Spectroscopy 

We obtained long slit low resolution spectra with FORS2 at 
VLT on Paranal on 2002 March 25th. The VLT uses an ac- 
tive optics platform to achieve high quality image. The see- 
ing conditions were excellent (~ 0'.'64) and we obtained high 
quality spectra. We used a 0'.'7 slit and the GR600I-I-25 and 
GR1200R-1-93 grisms, therefore covering a large wavelength 
range in the red part of the optical spectrum (Table The 
spectra were processed using a custom pipeline based on stan- 
dard procedures in IRAF. First, the two-dimensional images 
for the two separate CCDs which make up the FORS2 de- 
tector were independently overscan-, bias- and flat field- cor- 
rected. The lower chip image was then normalised to the me- 
dian value of the upper one, and the images for the two chips 
could then be merged using the fsmosaic tool of the FIMS soft- 
ware (FORS 1-1-2 FIMS manual, issue 2.6). The 1-D spectra 
were then extracted, and wavelength and flux calibrated in a 
standard manner, using the best calibration files provided by 
the VLT team. 

We obtained a HST-STIS low spectral resolution observa- 
tion, aiming for spatially resolved spectra of the two compo- 
nents. Observations occurred on 2003 June 26th. We used the 
G750L grating (0.525yum-1.300yum, 4.92A/pixel) with the U.'l 



In order to estimate the spectral type of DENIS-P J 185950.9- 
370632, we compare its optical spectrum to that of brown 
dwarfs from the field, and from the Upper Scorpius OB as- 
sociation (hereafter USco, iDelfosse et aljEoOSl iMartm et alJ 
2004). The age and distance of USco, resp ectively esti- 
mated at 5 Myr ([Preibisch & ZinneckeJ 19991) and 145 pc 
dde Zeeuw et al. l>1999'). are very close to those of R-CrA, and 
Ay is close to z ero, with local maxima reaching ~1.0 mag 
(see ICambresvl 1 999l) . Each reference spectrum was artificially 
reddened with different values of extinction (using the red- 
den task of IRAF) and then compared to DENIS-P J185950.9- 
370632 spectrum. Figures |6l and shows the best match ob- 
tained with the USco M8 brown dwarf DENIS-P J161916.5- 
234722.9, and the M8 field ultracool dwarf VB 10, both for 
an extinction of Ai/=0.5 mag. The correlation is very good 
and the two results agree perfectly. One can note that the K I 
and Na I doublets are stronger in VB 10, as expected be- 
cause of its higher gravity, while they are very similar in the 
case of DENIS-P J161916.5-234722.9. This is another hint that 
DENIS-P J 185950.9-370632 is young and is likely to belong 
to the R-CrA association, which age is very close to that of the 
Usco OB association. The next best match was with DENIS- 
P J 1556-2 106 (USco, M7) for a reddening of Av=1.0. We 
thus adopt a spectral type of M8+0.5 for DENIS-P J185950.9- 
370632 and an extinction of Av=0.5+0.3 mag. 

We applied the NICER (Near Infrared Colour Excess 
Revisited) technique ( Lombardi & Alves 2001) on the 2MASS 
catalog photometry to produce an extinction map of the 
R-CrA region, with a 2' resolution (Figure |8|i. DENIS- 
P J 185950.9-370632 lies next to a region of strong absorp- 
tion. The mean V band extinction in the 2' pixel that con- 
tains DENIS-P J185950.9-370632 is Ai/=3.7±0.4 mag ^. As 
explained above, the observed optical spectrum implies instead 
that DENIS-P J 185950.9-370632 is reddened by Ay -0.5 mag. 
This suggests that DENIS-P J 185950.9-370632 probably Ues 
somewhat on the near side of the cloud, or in a relative gap of 
the patchy extinction. 



' IRAF is distributed by National Optical Astronomy 
Observatories, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under contract with the National 
Science Foundation. 



- The 0.4 mag uncertainty quoted here is the measurement uncer- 
tainty on the mean extinction; the standard deviation of the extinction 
inside the 2' pixel is 1.8 mag 
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3.2. R-CrA membership 

3.2.1. Photometric distance 

After correcting for Ai/=0.5 mag, the magnitude and colour 
of DENIS-P J185950.9-370632 are / ~ 16.8 mag and 
I - J ma^ . The DUSTY theoretical isochrones 

dChabrier et a l.' 2000) associate to those parameters a distance 
of -150 pc at 5 Myr and 110 pc at 10 Myr. Corrected for bi- 
narity (assuming a flux ratio /2//1 =0.3 as measured in the 
F814W filter, see Tableland Section l33t . it corresponds to 
distances of ~170 pc and 125 pc. The corresponding masses 
are respectively -0.025 Mq and -0.035 M©. 

The late-M d\y a rfs sequence in Upper Scorpius 
jPelfosse et alJ l2003t iMartm et all l2004h provides a san- 
ity check: M7.5 dwarfs in USco have I between 16.0 and 
17.5 mag for M6.5 to M7.5 dwarfs, indicating for DENIS- 
P J 185950.9-370632 a photometric distance similar to that of 
Upper Scorpius, whose distance is close to that of the CrA 
complex. 

3.2.2. Proper motion 

As discussed in section lO DENIS-P J 185950.9-370632 has 
been detected at many epochs. Figure shows that HST/ACS 
position is clearly discrepant, which we suspect is due to 
underestimated uncertainties on the position obtained with 
HST/ACS, which lies over 3-cr away from any other Direct 
comparison with the HST/WFPC2 images confirms that the 
pointing is correct, and we suspect that the problem arises in 
the astrometric processing of the ACS data by the STSDAS 
pipeline. The DENIS, 2MASS, HST/WFPC2, HST/STIS, 
GSC 2.21 and USNO-Bl.O observations afl agree to within 
less than 1-cr for epochs spread over 18 years, confirming that 
the HST/ACS measurement is most suspicious. There are un- 
fortunately not enough other stars in the field of the HRC to 
perform meaningful and precise astrometric re-calibration. We 
derive an approximate motion for DENIS-P J 185950.9-370632 
by comparing the most precise of recent measurements, from 
2MASS, with the earliest position reported in the USNO-Bl.O 
(with similar uncertainties, -14 years earlier). This rough esti- 
mate of the proper motion, fia cos 5 - Qf.'Q yr ' + 0'.'013 yr ' 
and Us = -0'.'021 yr~ ' + 0^.^013 yr~\ agrees w ithin the (large) 
uncertainties with the lNeuhauser et all ( I2OOOI) value for the R- 
CrA region, fi^ cos 6 = 0'.'005 yr^^ and fig = -0'.'027 yr ' . This 
is consistent with membership of DENIS-P J 1 85950.9-370632 
in the star forming region, though the significance of the result 
is obviously not very high. 

3.3. Imaging: a close companion 

As shown in Figures El and DENIS-P J185950.9-370632 
is clearly elongated in the high resolution HST images, at 2 
epochs and in 6 different filters. We analysed the HST images 
for precise separations, position angles and flux ratios of the 
possible multiple system, us ing a custom-made program de- 
scribed in Bouy et al. (2003) and adapted here for ACS/HRC. 
Briefly, the PSF fitting routine builds a model binary using ten 



different PSF stars from several ACS/HRC images, and then 
perform a non-linear PSF fit of the observed image to deter- 
mine the best-fit values for the 3 free parame ters: separation, 
position angle and flux ratio. IBouv etal.lll2003 l) discuss the un- 
certainties and limitations of the algorithm in detail, but slight 
improvements since then have led to a much better understand- 
ing of the uncertainties and systematic errors t.Bouv..2004k , in 
prep.). 

Table |3l summarises the resulting best-fit binary parame- 
ters. For the three filters where the source is best resolved 
(F625W, F675W, F775W; see Figures |4| and H they are fully 
consistent, demonstrating that a binary star is an excellent 
model of the observations. The nominally discrepant param- 
eters in the other filters are from marginally resolved im- 
ages, either due to diffraction broadening at redder wavelength 
(WFPC2/PC F814W, ACS/HRC F850LP) or because of wider 
pixels (STIS/F28X50LP). They should therefore not be given 
much weight. Given the proper motion derived in section lT.2.21 
two independent objects should have moved apart by -0'.'042 
between the two observations. The lack of any apparent change 
is thus a strong indication that they form a common proper mo- 
tion pair The flux ratios indicate that the secondary is signifi- 
cantly fainter than the primary at all observed wavelengths. 

3.4. Spectral Analysis 
3.4.1. Spectral Features 

Tablel^lists the equivalent width (EW) of several spectral fea- 
tures measured in the different data sets. The only emission 
line in the HIRES spectrum is Ha. Other common emission 
lines such as He I at 667.8 nm, O I at 844.6 nm, and Ca 11 at 
866.2 nm are not detected, with upper limits on their equivalent 
width (EW) below 0.5 A. By direct integration of the line pro- 
file (using the splot IRAF task) we measure an Ho- equivalent 
width of 18+3 A in the HIRES specti'um, 18±3 A from the 
low-resolution STIS spectru m, and 17+2 A from the FORS2 
spectrum. According to Barr ado v Navascues & Martini ( l2003h 
EW(Ha)=47.5 A for an M7.5 spectral type implies that the 
Ha emission is caused by accretion, while weaker lines can 
be due to either accretion or chromospheric activity. The Ha 
line strength by itself is thus here insufficient to distinguish 
between an accretor and a chromospherically active star The 
Ha line however is relatively broad (Figure |5Jl, with a full 
width at 10% of peak intensity of 205±10 km s"'. Ha line 
widt hs above 200 km s~' in b rown dwarfs are due to accre- 
tion j-Tavawardhana et ai]l2003t) . and by this measure DENIS- 
P J185950.9-370632 is slighty above the limit and is likely 
an accretor If accreting, the mass accretion rate is probably 
rather low, given the modest Ha strength and the lack of op- 
tical veiling. Accretion rates below 10"^ Mq yr"' produce no 
measurable veiling (MuzeroUe et al. 2000) . The consistency of 
the EW measured at the three epochs suggests that we are mea- 
suring quiescent emission rather than variable activity. This is 
more consistent with steady accretion than with chromospheric 
activity, where strong and broad Ha lines are only observed 
during outbursts. Further observations should be made in order 
to confirm that DENIS-P J 1 85950.9-370632 is accreting or not. 
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In the region of the Li I resonance line (670.8 nm) the 
HIRES spectrum is rather noisy, but the line is well de- 
tected with EW(Li I)=0.41+0.08 A. This line also appears 
clearly in the FORS2 spectrum, as shown in Figure |6] with 
EW(Li I)=0.9+0.4 A. The resolution of the STIS spectrum is 
insufficient to isolate the Li I line. The K I line at 769.9 nm is 
also present in the HIRES and FORS2 spectra. As already men- 
tioned in section im Figure^shows that the K I and Na I lines 
of DENIS-P J 185950.9-370632 are much narrower and weaker 
than in the spectrum of a field M8 dwarf (VB 10), indicating a 
young age. 

3.4.2. Comparison with model spectra 

In order to estimate of the photospheric effective temper- 
ature, we compare the observed spectra with the DUSTY 
jAUard et al.ll200lh and NextGen dHauschildt etali.1999 ) at- 
mospheric models. The synthetic spectra were smoothed with 
a gaussian kernel matched to the slit size and resampled to the 
same grid as the observed spectra. The best fit model was de- 
termined by maximizing the cross-correlation with the dered- 
dened observed spectrum (a minimum adjustment, per- 
formed as a cross-check, gives identitical results). 

The free parameters for the fit are the effective temperature 
and the surface gravity (the metallicity was fixed to solar metal- 
licity; models for higher metallicities are not yet available, and 
lower metallicities gave significantly degraded agreement). The 
latest NextGen and DUSTY models give very similar surface 
gravities. For an age of 5 Myr and I-J~3.0, they respectively 
give log g ~3.8 and log g ~4.0, and for 10 Myr respectively 
log g ~4.l and log g ~4.0. The available DUSTY models have 
a 100 K grid step, and cover surface gravities log g ranging 
from 3.5 to 6.0 with a 0.5 interval. We only considered mod- 
els with effective temperature between 1500 K and 3000 K, 
and surface gravities between 3.5 and 4.5. The observed and 
synthetic spectra were both normalized to an integrated flux of 
unity prior to the analysis. 

Figure [TO] shows the best-fit DUSTY synthetic spec- 
trum, over-plotted on the observed spectra. For both sets 
of models (DUSTY and NextGen), the best fits are ob- 
tained for Te//=2600 K and log ^=3.5 (STIS spectrum), and 
Teff-27QQ K and log §=3.5 (FORS2 spectrum). These temper- 
atures correspond to that of the equivalent unresolved system. 
Attempts to account for the multiplicity by using two synthetic 
spectra at different effective temperatures did not produce a sig- 
nificantly improved fit. This is actually expected, since accord- 
ing to the DUSTY models the observed magnitude differences 
(from 0.60 mag in the HST F775W filter to 1.74 mag in the 
HST F625W filter) correspond to effective temperatures that 
differ by <200 K. This is close to the effective temperature res- 
olution of our adjustment, and a dual temperature fit is there- 
fore not warranted. The obtained effective temperature should 
on the other hand represent a reasonable estimate for both com- 
ponents. 

It is important to note that the model spectra do not match 
the observed ones very well over the present spectral range, 
perhaps because of their simplistic handling of the gravitational 



settling of dust. They probably underestimate the strength of 
molecular absorption bands like TiO, VO, but of also atomic 
lines like K I and Na I D. Since these features dominate the 
energy distribution in optical spectra of late M dwarfs, the pa- 
rameters that we obtain most likely suffer from systematic er- 
rors. They should on the other hand be much more reliable 
when used in a relative sense, and compared with other anal- 
yses based on the same atmospheric models. The low surface 
gravity of DENIS-P J185950.9-370632 (log g = 3.5) relative 
to field M8 dwarfs (log g~4.5), in particular, is a robust result. 

4. Discussion 

4.1. Infrared excess 

Figure ^2 compares the spectral energy distribution (SED) 
of DENIS-P J 185950.9-370632 to the DUSTY and NextGen 
models (for 5 Myr), as well as to observed field late-M dwarfs. 
The fluxes have been normalized to have an integrated lumi- 
nosity between and 1.65 yum (H band) equal to one. The 
choice of the H band as limit instead of, for example, the over- 
all SED, was made for the following reasons: first because the 
SED of late-M dwarfs peaks around this value (see the DUSTY, 
NextGen and field M dwarfs SED in the figure), and second be- 
cause in this wavelength range, the SED should be less affected 
by accretion-related continuum emission and infrared emission 
from the disk than at redder colours. 

DENIS-P J 185950.9-370632 shows a sti'ong excess over 
the models and the field dwarfs at wavelength greater than the 
H band. This wavelength range (A >1.6 //m) is on the Raleigh- 
Jeans tail of the spectrum, and models are expected to be very 
reliable. We did not find any published observations of field 
late-M dwarfs in the 5.0 jt/m-8.5 fim range, but the available 
photometric measurements in the R to L' bands already show 
that the field M dwarfs have a much lower flux between K and 
L than DENIS-P J185950.9-370632, independently demon- 
strating the infrared excess. This excess refers to the unresolved 
system, which is probably dominated by the primary (see sec- 
tion |33li. But since DENIS-P J 185950.9-370632 is a binary, 
we checked that its SED cannot be fitted by any meaningful 
combination of two synthetic SED, and the infrared excess is 
robust. 

4.2. An accreting close Binary 

As discussed in Section l331 DENIS-P J185950.9-370632 is 
very likely to be a binary, with a separation of ~0'.'060, 
or 7.8 A.U at its 130 pc distance. This separation is well 
within th e 0-20 A.U range observed for field brown dwarf 
systenis dBouv et alJl2()03HBurgasser et all 120031: IClose et alJ 
l2003t 'Gizis et al.' '2003V The statistical correction factor of 
1.26 from Fis cher & M arcv (1992), leads to a semi-major axis 
of a =9.8 A.U. As discussed in section |3 the photometry, 
the spectral distribution, and the optical low resolution spec- 
tra all indicate an age of 5~10 Myr for a distance of 130 pc, 
and an effective temperature of ~2 600 K. The small magni- 
tude difference (Table |3j indicates that the two components 
of the system must have fairly similar masses. According to 
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the DUSTY models, magnitude differences of -1.5 mag in the 
F675W filter and ~1.1 mag in the F814W filter correspond to a 
mass ratio of M^/M^ ~75%. Considering a total mass between 
0.020 M0 <M,o, <0.030 M© and an orbit with a semi-major 
axis of 9.8 A.U, this leads to an orbital period of ~ 170 yrs 



Aeronautics and Space Administration and the National Science 
Foundation. This work made use of data from the The Guide Star 
Catalogue II, which is a joint project of the Space Telescope Science 
Institute and the Osservatorio Astronomico di Torino. 



5. Conclusions 

The results presented in this paper demonstrate that DENIS- 
P J 185950.9-370632 is a young multiple system in the R-Cra 
star forming region. On the optical images, we find a clear elon- 
gation. After PSF subtraction on 3 epochs data we conclude 
that DENIS-P J 185950.9-370632 is very likely to be a common 
proper motion pair with a separation of ~0'.'060, close to the 
resolution limit of the instruments we used. The spectroscopy 
constrains the effective temperature to 2600~2700 K for a low 
surface gravity (log ^=3.5), consistent with a young age. This 
temperature corresponds to a total mass of ~0.030+0.010 Mq 
for an age ranging between 5 and 10 Myr. This mass is consis- 
tent with the photometry, and DENIS-P J 185950.9-370632 is 
therefore clearly substellar. Infrared excess observed with the 
ISOCAM in the 5.0-8.5//m band and the presence of a strong 
Ha emission as well as lithium absorption at 670.8 nm suggest 
a young age, the presence of circumstellar material, and a sub- 
stellar mass object. Added to a consistent preliminary estimate 
of the proper motion and consistent colours, these observations 
suggest that DENIS-P J 185950.9-370632 belongs to the R- 
CrA star forming region. From the magnitude difference be- 
tween the two components we estimate a mass ratio of ~75%. 
The estimated orbital period is about 170 yrs. The mid-infrared 
excess observed by ISO and the width of the Ho' emission sug- 
gest that it might be surrounded by a disk and slowly accreting. 
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Table 2. Observation log. 



Imaging 


Instnimeiit Filter 


Exp. 


Ualc VJUS. 


Pixel Scale 




Xime [s] 




L J 


WFPC2-PC F675W 


600 


19/09/2000 


a'045 


WFPC2-PC F814W 


300 


19/09/2000 


a'045 


ACS-HRC F625W 


1000 


24/09/2002 


a'027 


ACS-HRC F775W 


460 


24/09/2002 


a'027 


ACS-HRC F850LP 


340 


24/09/2002 


a' 027 


ISO-CAM 1 LW2 


4328 


20/04/1996 


6'.'0 


ISO-CAMl LW3 


4326 


20/04/1996 


6ro 


Spectroscopy 


Instrument Wavelength 


Exp. 


Date Obs. 


Dispersion 


Range [//m] 


Time [s] 




[A/pixel] 


VLT-FORS2 0.590-0.715 


600 


06/05/2002 


0.60 


VLT-FORS2 0.690-0.910 


200 


06/05/2002 


1.06 


HST-STIS 0.525-1.300 


4140 


16/06/2003 


4.92 


Keck-HIRES 0.667-0.895 


1200 


30/05/2000 


1.1 


Table 4. Spectral Features 






Feature Instrument 


A [nm] EW [A] 




Ha HST/STIS 


656.3 -18 


±3 




Ha VLT/F0RS2 


656.3 -17 ± 2 




Ha Keck/HIRES 


656.3 -18 ± 3 




Na I VLT/F0RS2 


818.3 3.7 + 0.3 




Li I VLT/FORS2 


670.8 0.9 


± 0.4 




Li I Keck/HIRES 


670.8 0.41± 0.08 
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Epoch 


R.A 


Dec. 


Uncert. 


Filter 


Mag. 


Source 


01/01/1985 


18 59 50.9 


-37 06 31.0 


±a'i 


Rl 


18.9 


USNO-Bl.O 










B2 


21.7 












R2 


19.9 












I 


17.3 




06/07/1989 


18 59 50.9 


-37 06 31.2 


±a'3 


F 


20.2 


GSC2.21 


20/04/1996 


18:59:50.7 


-37:06:28.0 


±6'.'0 


LW2 


11.6 


ISO 










LW3 


>10.7 




28/04/1999 


18:59:50.9 


-37:06:32.0 


±1'.'0 


I 


17.01 


DENIS 










J 


13.99 












K 


12.60 




25/06/1999 


18:59:50.9 


-37:06:31.3 


±a'i 


J 


13.98 


2MASS 










H 


13.10 












K 


12.56 




12/09/2000 


18:59:50.9 


-37:06:30.6 


±1'.'0 


F675W 


19.33 


HST/WFPC2 










F814W 


16.88 




24/09/2002 


18:59:50.8 


-37:06:27.6 


±1'.'0 


F625W 




HST/ACS 










F775W 














F850LP 






26/06/2003 


18:59:50.8 


-37:06:31.4 


±1'.'0 


Long pass 19.6 


HST/STIS 


Table 3. PSF fitting results 


Date Obs. 


Instr. 


Filter 


Sep.' 


["] 


PA' n 


AMag' 


23/06/2003 


STIS/CCD 


F28X50LP 


elong 


ated"^ 


elongated"^ 


elongated" 


24/09/2002 


ACS/HRC 


F625W 


a.'065+a'OOl 


279.2±0.1 


1.74+0.06 


24/09/2002 


ACS/HRC 


F775W 


ff.'057±a'0005 


279.1 ±0.1 


0.66±0.05 


24/09/2002 


ACS/HRC 


F850LP 


elong 


ated 


elongated 


elongated 


12/09/2000 


WFPC2/PC 


F675W 


a'066±0'.'003 


283.8±1.2 


1.30±0.11 


12/09/2000 


WFPC2/PC 


F814W 


a'059+a'003 


271.8+1.2 


1.10+0.11 



' The uncertainties reported here are 1-cr uncertainties as explained in iBouv et alj(i2003i) . 

^ As explained in the text, the PSF fitting program did not give good enough results, although the object is clearly elongated. 

Young, E. T. 1992, ApJ, 397, 520 
Wilking, B. A., McCaughrean, M. J., Burton, M. G., et al. 1997, 

AJ, 1 14, 2029 
York, D. G. e. a. 2000, AJ, 129, 1579 
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Fig. 1. The environment of DENIS-P J 185950.9-370632 in the CrA molecular cloud complex. The cross marks the position of 
the DENIS objects on this DSSl/STSci J plate. The field of view is 50'x 35'. North is up and East is left. 
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Fig. 2. Finding charts of DENIS-P J185950.9-370632in the I and J filters (DENIS images). The field of view is 3.5'x 3.5', North 
is up and East is left. 
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Fig. 3. HST-WFPC2/PC image obtained in the F675W and F814W filters (composite image) on 2000 September 12th. The over- 
plotted circles represent the position of DENIS-P Jl 85950.9-370632 at different epochs and with different instruments. The radius 
of each circle corresponds to the uncertainty of that measurement. Although the different measurements are not fully consistent, 
no other source than DENIS-P Jl 85950.9-370632 can be associated in the image. 
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Fig. 4. Contour plots of the ACS/HRC (F625W) and WFPC2/PC (F814W) images. The WFPC2/PC and ACS/HRC image have 
identical scales for easier comparison. The object is clearly elongated on both images, but is better resolved on the higher 
resolution and better sampled ACS image. 
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Fig. 5. Surface plots showing the results of the non-Unear PSF fitting on the ACS/HRC (F625W filter) and WFPC2/PC (F814W) 

images. Amplitudes are normalised. The sky background has been subtracted. The figure shows the images obtained with 
HST/HRC and WFPC2/PC, one of the PSF stars, the companion appearing after PSF subtraction, and the residuals after sub- 
tracting the modelled binary system. The pixel scale of the ACS/HRC is almost twice finer than that of the WFPC2/PC. 
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Fig. 6. Comparison of DENIS-P J185950. 9-370632 FORS2 spectrum (smoothed via a boxcar with a width of 5 pixels) with 
the spectrum of the Upper Scorpius member M8 dwarf DENIS-P J161916. 5-234722. 9 artificially reddened with Av=0.5 mag. 
The match is very good. The inset box shows a zoom of DENIS-P J 185950.9-370632 spectrum around the Li 1 absorption. 
Spectrum of DENIS-P J161916.5-234722.9 from iMartm etaD (|2004 . Fluxes have been normalized at the pseudocontinuum 
level at 840.0 nm. 
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Fig. 7. Comparison of DENIS-P J185950.9-370632 FORS2 spectrum (smoothed via a boxcar with a width of 5 pixels) with 
the spectrum of the field M8 dwarf VB 10, artificially reddened with Av=0.5 mag. The match is very good. The inset boxes 
show zooms of the two spectra around the K I and Na I doublet. It shows clearly that these two doublets are stronger in VB 10, 
as expected because of its higher gravity. Spectrum of VB 10 from .Martin et aL (.1999.) . Fluxes have been normahzed at the 
pseudocontinuum level at 840.0 nm. 
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Fig. 8. Extinction map of the R-CrA Region around DENIS-P J185950. 9-370632 (indicated by white cross). Isocontours of 
Ay=0.0, 1.0, 2.0, 3.0 and 4.0 mag indicate the scale. 
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Fig. 9. HIRES spectra of DENIS 1859-37 (solid line) and Gl 406 (dotted line) in the H„ region. Gl 406 is a typical chromo- 
spherically active M6 dwarf. DENIS 1859-37 has a much broader and asymmetric emission line, indicating that it is probably an 
accretor. 
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Fig, 10. STIS and FORS2 spectra compared with models. Top panel: STIS spectrum of DENIS-P J185950.9-370632 compared 
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Fig. 11. Spectral energy distribution of DENIS-P J 185950.9-370632 compared with DUSTY and NextGen models (5 Myr) and 



